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Possibility of Solid-Fluid Transition in Moving Periodic Systems 
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The steady sliding state of periodic structures such as charge density waves and flux hue 
lattices is numerically studied based on the three dimensional driven random-field XY model. 
We focus on the dynamical phase transition between plastic flow and moving solid phases 
controlled by the magnitude of the driving force. By analyzing the connectivity of comoving 
clusters, we find that they percolate the system under driving forces larger than a certain critical 
force within a finite observation time. The critical force, however, logarithmically diverges with 
the observation time, i.e., the moving solid phase exists only within a certain finite time, which 
exponentially grows with driving force. 
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The collective transport phenomena of condensed mat- 
ter with random pinning attract much attention from 
the viewpoints of solid state physics, nonlinear dynamics 
and statistical mechanics. There are numerous systems 
which belong to this class of dynamics, e.g., charge den- 
sity waves (CDWs),^ flux line lattices (FLLs),^ colloidal 
lattices'^ and Wigner crystals.^' ^ These systems have spa- 
tial periodicity in the absence of random pinning, which 
modifies the periodic order and pins the system. Under 
a driving force larger than a certain threshold value, the 
system starts moving and shows highly nonlinear conduc- 
tion. This is the depinning transition, which has been in- 
vestigated extensively. Recently another attractive topic 
of these systems has been the dynamical phase transition 
between two nonequilibrium steady states above the de- 
pinning threshold field. In the "ordered phase" the 
local DC velocity is uniform and the spatial periodicity 
has a quasi long range order. Such a state in FLL sys- 
tems is called moving Bragg glass. ^'^ In the "disordered 
phase" , which is often called the "plastic flow phase" , the 
motion is spatially nonuniform and the periodic order is 
destroyed. This transition is considered to be induced by 
a change in driving force or pinning strength. Although 
there are a lot of experimental studies^*^' and numeri- 
cal simulations^^"^^ that show an evidence of a dynamical 
melting transition between moving Bragg glass and fluid 
phases, the existence of such an ordered phase is still an 
unanswered question. 

In such discussions, it is assumed implicitly that two 
types of order, namely, the spatial periodic order and 
the local DC velocity order, are established simultane- 
ously. To date, the former has been mainly discussed. 
These two orders, however, are independent in principle. 
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e.g., the "moving glass state", in which the spatial peri- 
odicity is destroyed but frozen in time so the local DC 
velocity is uniform, is possible. In this article we focus on 
the dynamical property, i.e., the uniformity of the local 
DC velocity. It is closely related to the plastic deforma- 
tion because at the boundary between the domains that 
have different velocities, the local strain increases with 
time and tearing occurs. We discuss the phase transition 
between the "plastic flow phase" and the "moving solid 
phase" , which are distinguished by the existence of the 
local DC velocity order. 

We perform numerical simulations based on the three 
dimensional driven random-field XY model. ^^"^"^ It is a 
modified version of the intensively investigated elastic 
manifold model, such as the Fukuyama-Lee-Rice model 
for CDWs.^^'^^ The density of a periodic structure is ex- 
pressed as p{r, t) ~ pq cos (q • [r — u(r, t)]) . Here, q's 
are fundamental reciprocal lattice vectors. Higher har- 
monics are ignored here. u(r, t) is a deformation field and 
the phase field 0q(r, i) = q • u(r, t) is often employed as a 
degree of freedom. In this article, only the case of a single 
phase field, which denotes the driving direction compo- 
nent of deformation, is treated. It is sufficient in the case 
that periodicity is one dimensional as CDWs in NbSes. 
We consider, however, that the essential feature of the 
dynamics of higher dimensional periodic systems such as 
FLLs would be captured. The periodicity of the structure 
is related to the phase coherence and experimentally ob- 
servable transport quantities, such as electric current for 
CDWs and voltage drop for FLLs, are proportional to 
the phase velocity 9{y, t). 

The elastic manifold model treats internal interaction 
using the elastic energy / dr{\/9)'^, which becomes a har- 
monic coupling, j) (^i ~ ^j)^/^, in a lattice model. In 
order to treat plastic deformation, we replace this har- 
monic coupling with a sinusoidal one, 1 — cos{9i — 9j). 
They are equivalent in the limit where phase differ- 
ences become zero. Here, the indices i's denote semi- 
macroscopic domains fixed in the space in which phase 
coherence is always held. This sinusoidal coupling in- 
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duces maximum restoring force, i.e., yield stress and al- 
lows plastic deformation, so-called phase slip. Phase slip 
is a process in which the phase difference between neigh- 
boring domains increases or decreases by 27r. It results 
in no change in coupling energy. 

The overdamped equations of motion for the phases of 
domains, 6i's, are as follows. 

Oi = --Y^ sm{0^ - 0j) - sin(0, - P^) + f (1) 
j 

Wc choose the units that both of a pinning strength and 
a dissipation coefficient equal unity. The first term on the 
right hand side refers to the interaction with neighboring 
z domains. The second term denotes the a random pin- 
ning force and /3j's are given as uniform random numbers 
between and 27r. / is a uniform driving force. 

Strogatz et al. analyzed this model by the mean field 
approximation and found a discontinuous transition by 
changing the external field. There are three regimes, a 
pinned static state and a homogeneously moving state, 
and a bistable regime between the two. Huse performed 
numerical simulations of this model in three dimen- 
sions.^^ He investigated the Lyapunov exponent and ve- 
locity coherence and found a transition between spa- 
tially uniform temporally regular motion and nonuniform 
chaotic motion by changing pinning strength. These mo- 
tions are related to moving solid and plastic flow, respec- 
tively. In this article, we analyze the / dependence of the 
dynamics of this model systematically paying special at- 
tention to the dependences on system size and observa- 
tion time. 

We numerically solve eq. (1) by the fourth-order 
Runge-Kutta method. The discretized time step is set at 
2-!t/8{J + f). The domains are put regularly on the sim- 
ple cubic lattice in three dimensions. We call this unit 
"site" instead of "domain" hereafter. A periodic bound- 
ary condition is imposed. There are two independent pa- 
rameters, coupling constant J and driving force /. In 
this article we show mainly the results for J = 1.0. All 
phases are set at the unique value in the initial state. 
Physical quantities are calculated after some precursor 
running (typical time is 12900) for the relaxation to the 
steady state. An initial state is sometimes substituted by 
the final state for the simulation with a slightly larger /. 
Thermal fluctuation is not taken into account. Simula- 
tions are performed with some samples that have differ- 
ent set of /3i'ii. The numbers of samples whose linear sizes 
are 16,32,64 and 128 are 32,16,8 and 4, respectively. 

In the left panel of Fig. 1 a spatial configuration of the 
local DC velocity, ij^q = {9i(t))T, is shown. Here, {■ ■ ■)t 
denotes time averaging for the observation time T and 
the resolution of DC velocity is given by 2it/T. Note that 
we show the results of two dimensional systems, which 
show qualitatively similar behaviors to three dimensional 
systems but show the domain structure of u;^^ more 
clearly. In the right panel we show the difference in uj^^q, 
l^'^Dcl ~ I'^DC ~ '^Dcl' ^'-^^ each bond between neigh- 
boring sites. When lAwp-^j > 2'jt/T, phase slip occurs on 
the bond at least once during the observation time T. 
We focus on this phase slip process to discuss the spatial 




Fig. 1. Spatial configuration of ludc on site (left) and Aludc on 
bond (right) in a two dimensional sample. uj-ocT/2tt is plotted 
with gray scale from 2730 to 2820 as color changes from white 
to black and log(|Au)Dc|r/27r) changes from logO.25 to log 1000. 
Connected bonds are plotted in white. Result of two dimensional 
system with 256^ sites, T « 25700, J=1.0 and / = 1.0 (/c(T) = 
1.06). 

correlation of motion instead of the direct spatial corre- 
lation of DC local velocity. -'^■^ We define such a bond as 
a "disconnected" bond. Otherwise, if |A(x;dc| < 2it/T, 
the bond is "connected" and two sites belong to the 
same cluster. When both of the pair sites are pinned, 
i.e., < 27r/T, IAwdcI is less than 2it/T but 

we regard such a bond as disconnected, which forms 
nonmoving solid clusters. We analyze a bond percola- 
tion transition by controlling driving force. The moving 
solid phase is characterized by an infinitely large cluster, 
which is made from connected sites without phase slip. 
The percolating phase is thus the moving solid phase. 
The driving force in Fig. 1 is slightly below the critical 
point and a fractal domain structure appears. 

In order to perform finite size scaling, we divide the 
system into subsystems, whose linear dimension L is 
smaller than that of a real sample Lmax- Then we deter- 
mine that percolation occurs if a certain cluster reaches 
the two opposite sides of each subsystem. The statistics 
of subsystems and samples yields percolation probabil- 
ity P{f,L,T), which monotonically grows with / and 
becomes smaller as T increases. 

The reason wc perform finite size scaling in terms 
of L and not Lmax as usual is the following. We have 
found that the whole system with a finite imax even- 
tually falls into a periodic motion above some thresh- 
old driving force. Therefore, for a given driving force, 
we choose to work on a system with sufficiently large 
Lrnax to eliminate this artificial periodic motion and an- 
alyze the connectivity of clusters at various length scales 

-£'(< i'max)- 

In Fig. 2, P(/, L, T)'s for various L's and T's are plot- 
ted as functions of /. We consider the large L limit 
at fixed T first and the T dependence next. Percola- 
tion probability grows with / and its shape comes closer 
to that of a step function of / as L increases. Finite 
size scaling can be performed in the same way as the 
stochastic percolation. The curves for different L's con- 
verge on the universal one as the driving force is scaled 
as / = [(/ — fc)/fc]L^^'' with the suitable critical force 
fc and the critical exponent ly. The correlation length di- 
verges as |(/ — fc)/fc\~'^- A good conversion is obtained 
for each T as shown in Fig. 2 then /c(r)'s and i^(T')'s are 
obtained. The fraction of connected bonds at the critical 
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Fig. 2. Raw data of percolation probabilities for various sizes and 
observation times (top) and result of finite size scaling (bottom) . 
The X-axis of the latter is shifted by a constant which depends on 
T and each origin is indicated by a dotted vertical line. Result 
of three dimensional samples with Lmax = 64 and maximum 
observation time T ?s 3200. 



point is 0.10, which is smaller than that for the stochastic 
percolation, 0.2488,^^ due to the attractive correlation of 
connected bonds. 

The magnitude of the critical driving force for the per- 
colation transition obtained from the above analysis de- 
pends on T. The T dependence of /c for several Umax's 
is shown in Fig. 3. It is expressed as 

/,(r) = /oin(rAo), (2) 

for a long T. Here, /o and to are constants. We expect 
that eq. (2), which is consistent with the idea of shaking 
temperature*" as mentioned later, holds for a sufficiently 
long T and large Lmax and the reason for the deviation 
between the present data and cq. (2) for a smaller T or 
a smaller Lmax is considered as follows. For T < 2000, 
the fc's obtained from the simulations are a little larger 
than that expected from cq. (2). This discrepancy can be 
due to the existence of pinned or very low velocity sites, 
which cannot be candidates of the percolating cluster. 
They arc defects of the percolation transition, of which 
effects are not taken into account in eq. (2). The percola- 
tion transition is not so sensitive to the existence of such 
defects when they are rare and isolated. They are not a 
minority, however, below / = 0.80, e.g, the fractions of 
such sites with wdc < 0.050 are 0.51, 0.090 and 0.00054 
for /=0.75, 0.80 and 0.85, respectively. (Here, the local 
depinning threshold force^^'^° equals 0.10.) The critical 
fraction of the connected bonds becomes larger in the 
lattice with such defects and fc becomes higher. These 
defects decrease rapidly with / and have less effect on a 
higher for a longer T. 

Another disagreement occurs at a longer T, which is 
caused by the finiteness of the real sample size Lmax- 
For a finite imax, fc exhibits saturating behavior to a fi- 
nite value above a certain T, which increases with Lmax- 
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Fig. 3. f — T phase diagram. Phase boundary between plastic 
flow and moving glass is drawn by / = fc{T). Note that the 
horizontal axis is in a logarithmic scale. The dotted line indicates 
/ = 0.0251nt -1-0.59. 
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Fig. 4. Observation time dependence of critical exponent. 



As mentioned before, this is caused by the falling of the 
system into limit cycle motion when the phase coher- 
ence length becomes comparable to .Lmax- For T < 51500 
{fc < 0.870), the results for L = 64 and 128 hardly differ 
and they are considered to show fdT) for imax = oo. 

From these discussions eq. (2) is expected to hold up 
to the infinite T. This means that fc{T) diverges loga- 
rithmically as T approaches infinity and the moving solid 
phase does not exist in the long time limit based on the 
present definition. The supplemental simulation indicates 
that this behavior does not change for a system with a 
stronger coupling. 

From another point of view, eq. (2) is regarded as a 
type of "phase boundary" between the plastic flow and 
moving solid phases in the force-time plain (See Fig. 3). 
Considering the observation with fixed /, the crossover 
time r(/) is obtained as 



T(/)=toexp(///o). 



(3) 



The system behaves as if it were a moving solid in an 

observation time shorter than r(/). Beyond this time, 
cracks of plastic deformation, which are sheets of phase 
slip bonds, propagate to the macroscopic scale and fluid 
like property is revealed. 

Equation (3) can be regarded as a thermal activation 
process, r oc exp{V/kBTcs), if an effective temperature 
proportional to /~^(~ "-^dc) supposed. The inverse of 
/o is then proportional to the potential barrier V. This is 
consistent with the idea of "shaking temperature" pro- 
posed by Koshelev and Vinokur.'"' 

Next, we discuss the universality of the percolation 
transition for different observation times. In Fig. 4, the T 
dependence of the inverse of the critical exponent u is 
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Fig. 5. Relationship between maximum cluster size and linear di- 
mension of subsystem size for several observation times. 



shown. There is a tendency for large -Lmax's, 64 and 128, 
where has T-independent value « 0.8. The inverse 
of V becomes smaller and approaches zero as T increases. 

= means that P{f,L,T) does not depend on L, 
which manifests that the system falls into the limit cy- 
cle motion, i.e., a finite size effect. Such effect appears 
more clearly in the critical exponent than in the critical 
force, because the former is determined by the behavior 
in the whole critical regime. So we expect the percolation 
transition discussed here is universal with respect to the 
observation time. 

This imivcrsality is confirmed by the property of the 
percolating cluster at the critical point. Cluster size is 
identified with the number of contained sites and we de- 
fine Sc{L, T) as the size of the maximum cluster in a sub- 
system at / = fc{T). In Fig. 5, Sc(L,T)'s for several T's 
are plotted as functions of subsystem size L. They are ex- 
pressed as Sc{L, T) K, O.SSL^'^^ with a fractal dimension 
of 2.45 and show little T (or fciT) ) dependence. The 
large deviation of the data for Lniax=64 at the longest T 
is due to starting the limit cycle motion. 

The universality of the transition for different T's 
means that the fluctuation of DC velocity has a scale- 
less spatial pattern, which depends on neither / nor 
the phase coherence length, which grows with /, if one 
chooses proper time scale, i.e., sees the DC velocity in 
the resolution of 27r/T(/). 

In conclusion, we numerically investigate the possibil- 
ity of the dynamical phase transition between plastic flow 
and moving solid phases, which are distinguished by the 
existence of the long range order of local DC velocity. By 
analyzing the percolation of no-phase slip bonds and its 
observation time dependence, we found that the moving 
solid phase becomes unstable in a finite lifetime. 

The condition for the connected bond, that no phase 
slip occurs eternally, may seem too strict. For example all 
bonds necessarily take phase slips if thermal fiuctuation 
exists. It is important, however, to note that the local 
symmetry limy^oo ^^^dc = is destroyed from the be- 
ginning due to the random field. The bond steadily takes 
a phase slip in the same direction no matter how rarely 
it occurs. 

The crossover time t(/) can be defined clearly; the sys- 
tem behaves like a moving solid in a shorter time scale 
than t(/) and macroscopic plastic deformation occurs 
beyond t(/). This crossover time increases exponentially 
with the driving force. Its characteristic scale of growth 



/o = 0.025 is very small compared with other scales such 
as pinning strength (=1), therefore the crossover time in- 
creases very rapidly in the narrow region of / and over- 
comes the macroscopic time scale. This is a possible rea- 
son the moving solid phase is observed in experiments. 
The situation is similar to the case of structure glasses, 
whose viscosity grows quite large and they show slow 
dynamics, then it is hard to distinguish whether an equi- 
librium phase transition exists. 

We focused on the macroscopic plastic deformation 
and distinguished between plastic fiow and moving solid. 
This stance is different from the conventional interest 
in the liquid-crystal(Bragg glass) transition. Although it 
seems natural that these transitions occurs at the same 
time, the absence of the moving solid phase discussed 
here is not immediately related to the absence of the 
long range periodic order. For example the spatial phase 
order in a long span is possible in the plastic flow phase if 
the propagation of plastic deformation along the domain 
boundary is temporally localized and leaves no change 
before and after it. On the other hand we see that the 
saturation of phase correlation length to the system size 
results limit cycle motion and plastic deformation is sup- 
pressed. Then if the transition from liquid to crystal 
or Bragg glass, which is not observed in the range of 
our simulation, occurs at finite /, the present transition 
would happen at the same time. 

The numerical calculation was performed on the Hi- 
tachi SR8000 at the supercomputer center, ISSP, Uni- 
versity of Tokyo and the present study is financially 
supported by a Grant-in-Aid for Scientific Research 
(15540370) from the Japan Society for the Promotion 
of Science. 
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